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We present results from Monte Carlo simulations of the one-dimensional Ising spin glass with 
power-law interactions at low temperature, using the parallel tempering Monte Carlo method. For 
a set of parameters where the long-range part of the interaction is relevant, we find evidence for 
large-scale dropletlike excitations with an energy that is independent of system size, consistent with 
replica symmetry breaking. We also perform zero-temperature defect energy calculations for a range 
of parameters and find a stifi'ness exponent for domain walls in reasonable but by no means perfect 
agreement with analytic predictions. 
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I. INTRODUCTION 

Two main theories of the spin-glass state have been 
proposed: the "droplet picture, "^2*^4 and the replica 
symmetry-breaking (RSB) picturei^iSiL^ The droplet pic- 
ture states that the minimum-energy excitation of linear 
dimension I containing a given spin (a droplet) has an 
energy proportional to l^, where 9 is positive (assuming 
that the transition temperature Tc is finite). It follows 
that in the thermodynamic limit, excitations that flip a 
finite fraction of the spins cost an infinite energy. In ad- 
dition, the droplet picture states that these excitations 
are fractal with a fractal dimension ds < d, where d is the 
space dimension. By contrast, RSB follows Parisi's ex- 
act solution of the Sherrington-Kirkpatrick (SK) infinite- 
range model and predicts that excitations involving a fi- 
nite fraction of the spins cost a finite energy in the ther- 
modynamic limit, and are space filling^ i.e., dg = d. 

There have been several numerical 
attempta^i^Qi^i^^i^i^'^i^^i^^ to better understand the 
nature of the spin-glass state for short-range spin 
glasses. The data, which are on small system sizes, are 
consistent with a picture in which the "stiffness expo- 
nent" for dropletlike excitations is different from that 
for domain-wall excitations, although it is still unclear 
if the difference is due to the limited range of system 
sizeSfii or if it persists in the thermodynamic limit. 
Recently Krzakala and MartinjiS as well as Palassini 
and Youngfii find evidence for an intermediate picture 
(called "TNT" which stands for trivial-nontrivial) in 
which only a finite amount of energy is needed to 
generate system-size excitations in the thermodynamic 
limit, but their surface has a nontrivial fractal dimension 
less than d. 

In addition to finite-T Monte Carlo simulations, 
there has been a substantial amount of numerical 
work studying the energy of domain walls at T = 
0,18,19,20^21,22^23,24^25 j,^^ energy of the domain wall is 



found to vary as L where L is the system size and 
9 is positive (for systems with Tc > 0). The droplet 
theory assumes that the stiffness exponent for domain 
walls ^domain-wall is the Same as the stiffness expo- 
nent for dropletlike excitations ^droplet- In addition, 
^'domain-wall > 0. This is in contrast to the RSB and 
TNT pictures where the exponents are expected to be 
different and ^droplet = 0. To simplify the notation(2& 
from now on we will denote ^domain-waii by 9 and 6'dropiet 
hy9'. 

In this work we consider the question of whether there 
are one or two stiffness exponents for the case of Ising 
spin-glass models in one-dimension in which the interac- 
tions fall off with a power a of the distance. Depending 
on (7, the system can have either a finite transition tem- 
perature Tc (with the long-range part of the interaction 
relevant), or Tc = 0, in which case it can be either in 
the short-range or long-range universality classi^iSi The 
advantages of this system are the following: (i) a large 
range of system sizes L can be studied, (ii) there are an- 
alytical predictions for the stiffness exponent within the 
droplet picture, and (iii) within a single model, simply 
by changing a parameter, one can cover the full range of 
behavior, from the absence of a spin-glass phase to the 
infinite-range SK model. We perform both Monte Carlo 
simulations at low T (which indirectly gives 9'), for a 
value of a that has a finite Tc, and ground-state calcula- 
tions for a range of values of a corresponding to all three 
regimes: Tc > 0, Tc = (long range), and Tc = (short 
range). We use the parallel tempering22i2S approach for 
both the finite-T and T — studies. 

In Sees, m mil and IIVI we introduce the model, ob- 
servables, and details of the Monte Carlo technique, re- 
spectively. In addition we describe the phase diagram 
of the model in detail. The critical region of the disor- 
dered long-range Ising spin glass is analyzed in Scc.lVlin 
order to determine the critical temperature of the sys- 
tem and to ensure that in Sec. IVII where we study the 
system at finite T, the temperature is sufficiently far be- 
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FIG. 1: Geometry of the one-dimensional Ising chain: the 
Ising spins are placed equidistantly on a ring in order to en- 
sure periodic boundary conditions. Their distance Tij is de- 
termined by calculating the geometric distance between the 
different sites of the circle of radius R and circumference L 
(not all bonds are shown for clarity). 



low Tc such that critical-point fluctuations do not affect 
the data. The zero-temperature results are described in 
Sec. lVIII and our conclusions are summarized in Sec. lVIIll 



II. MODEL 

The Hamiltonian for the one-dimensional long-range 
Ising spin glass with power-law interactions is given by 




T~L — ^ ^ JijSiSj , 



(1) 



FIG. 2: Sketch of the phase diagram in the d-a plane for 
the spin-glass state of the disordered long-range Ising model 
with power-law interactions following Ref.|3- The light shaded 
region is where there is both a finite Tc and the spin-glass 
state is controlled by the long-range part of the interaction 
(i.e., ^LR > dsn and Slr > 0). The thick solid line separates 
the region of short-range behavior from that of long-range 
behavior, i.e., it corresponds to 6sr = ^lr, and is denoted 
by o-c(d), see Eq. (|7J. The thick dashed line separates regions 
where Tc = from regions where Tc > 0, i.e., it corresponds to 
^ = 0, where 6 refers to the greater of ^sr and 0lr. The dark 
shaded region is where there is no thermodynamic limit unless 
the interactions are scaled appropriately by system size. The 
calculations are performed for d = 1 (marked by a horizontal 
dotted line), for which (i) add) = 2, (ii) S > for cr < 1, 
and (iii) the thermodynamic limit does not exist for ct < 1/2 
(unless the interactions are scaled by a power of the system 
size). These values of a are marked. The simulations at finite 
T are performed for a — 3/4, d = 1, which is indicated by a 
cross. 



where the sites i lie on a ring of length L, as shown in 
Fig. n and 5*^ = ±1 represent Ising spins. We place 
the spins on a ring in order to ensure periodic boundary 
conditions. The sum is over all spins on the chain and 
the couplings Jij are given by 



(2) 



where the are chosen according to a Gaussian distri- 
bution with zero mean and standard deviation unity: 



1 



cxp(-ey2) . 



(3) 



The constant c{a) in Eq. Q is chosen to give a mean-field 
transition temperature = 1, where 



(4) 



Here 



denotes an average over disorder. 



The distance between two spins on the chain r.y is 
determined by = 2i?sin(a/2), where R is the radius 
of the chain and a is the angle between the two sites on 
the circle (see Fig. QJ. The previous expression can be 
rewritten in terms of L and the positions of the spins to 
obtain 
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The long-range Ising spin glass with power-law inter- 
actions has a very rich phase diagram in the d-a plane. 
This is summarized in Fig. |21 which is based on the work 
by Bray et alm^ and by Fisher and Husa^ who present 
a detailed analysis of the role of long-range interactions 
within the droplet model. Spin-glass behavior is con- 
trolled by the long-range part of the interaction if cr is 
sufficiently small, and by the short-range part if a is suffi- 
ciently large. More precisely, one has long-range behavior 
if the stiffness exponent of the long-range (LR) univer- 
sality class, is greater than that of the short-range 




FIG. 3: Stiffness exponent ^ as a function of a for d > 
di, where di is the lower critical dimension for short-range 
interactions. The dark shaded region corresponds to a < d/2 
where Tc diverges unless the interactions are scaled by a power 
of the system size. The light shaded region is where both 
Tc > and the behavior is governed by the long-range part 
of the interaction. 
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FIG. 4: Same as Fig. El but for d < dj. 



the zero-temperature calculations we use a range of a 
between 0.10 and 3.00. From Figs.|21and^ we see that 
this includes all the possible types of behavior: (i) Tc > 0, 
(ii) Tc = (LR), (iii) Tc = (SR), and (iv) SK-like 
behavior for a 0. 



(SR) universality class, 6'sr, and vice versa. In addition, 
there is an exact result for 6'lr,, namely^^iSi 

^'lr = d - cr, (6) 

so long-range behavior occurs if 

a<ac.{d)^d^esR{d) . (7) 

Equation @ indicates that critical exponents depend 
continuously on a in the long-range region, even though 
they are independent of a in the region controlled by the 
short-range part of the interaction. 

The condition for a finite-temperature transition is 
> 0, where 9 refers here to the greater of 0SR and 0lr- 
For the short-range model, there is a finite-temperature 
transition (i.e., 6'sr > 0) for d larger than the lower crit- 
ical dimension d/, which is found numerically to lie be- 
tween 2 and 3,i8ii9,20^30,3i Figs. 01 and H we show the 
expected variation of 9 with a both for d > di and d < di. 

For a < d/2 the model does not have a thermodynamic 
limit (Tc diverges) unless the interactions are scaled with 
an appropriate (inverse) power of L, i.e., c{a) — > for 
L ^ oo. Note that (7 = corresponds to the SK model. 

For d—1, the case studied in our numerics, one has^^ 

^sr(I) = -1, (8) 

and so from Eq. (Q, CTc = 2, in agreement with rigorous 
results by van Enter and van Hemmen.^^ This situation 
corresponds to Fig. ^ with ^SR = ^1- To obtain a finite 
Tc m d — \, one needs 1/2 < cr < 1 and, in the finite- 
T simulations, we choose a — 0.75 (see Fig. (J}. For 



III. OBSERVABLES 

In our finite-T simulations we primarily focus our at- 
tention on the spin overlap defined by 

i=\^jZSfSl (9) 

1=1 

where "a" and "/?" refer to two copies (replicas) of the 
system with the same disorder. In order to determine the 
critical temperature of the system, we study the Binder 
ratic'^^ g defined by 

Here, (• • •) represents a thermal average, v is the correla- 
tion length exponent, and Tc is the critical temperature 
of the system. Because g is dimensionless, curves for 
g{T,L) cross at the same Tc for all L. We also compute 
the spin-glass susceptibility Xsg ' where 

Xso - L[{q')U = L^-^'X[L""{T Tc)] , (11) 

and rj describes the power-law falloff of the spin-glass 
correlations at Tc. 

IV. EQUILIBRATION 

For the finite-T simulations, we use the parallel tem- 
pering Monte Carlo method^^'^^ as it allows us to study 
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FIG. 5; Average spin overlap [(g)]av as a function of Monte 
Carlo sweeps Ns-weep, which each of the replicas perform, av- 
eraged over the last half of the sweeps for different system 
sizes L. Note that the data appear to be independent of the 
number of sweeps. The data shown are for T = 0.10, the 
lowest temperature studied. 



larger systems at lower temperatures than is possible 
with conventional "single-spin-flip" Monte Carlo meth- 
ods. We test equilibration by the traditional technique 
of requiring that different observables are independent 
of the number of Monte Carlo steps iVgwocp- Figure [S] 
shows data for the spin overlap as a function of Monte 
Carlo steps for different system sizes. We equilibrate by 
repeatedly doubling the number of Monte Carlo sweeps 
until the last three agree within error bars. 

We also checked the equilibration of some of our 
data using a generalization of a test that was developed 
earlier— for nearest-neighbor models. For a Gaussian 
distribution of bonds it is straightforward to show, by 
integrating by parts with respect to the bonds, that the 
energy per spin, U = -N~'^Y.{i,j)[Jtj{SiSj)]av, can be 
related to a generalized "link overlap" qi , defined by 



Qi 



2 ■^-^ [Jij]a,v 



as follows: 



2T\U\ 

(TMF)2 



(12) 



(13) 



Note that this definition of qi is the natural generaliza- 
tion, to long-range interactions, of the link overlap dis- 
cussed in earlier workiiti^ on nearest-neighbor models, 
and that qi is normalized so that it tends to unity for 
T ^ where {S.Sjf 1. 

As discussed in Ref. 15, the two sides of Eq. H13|) are 
only equal to each other in equilibrium and approach the 
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FIG. 6: Test for equilibration using Eq. as discussed in 
Ref. 113 and the text. The data are for L = 128, T = 0.1. 
The figure shows both the link overlap qi defined in Eq. I12II 
and qi (energy) which is the RHS of Eq. , for an increasing 
number of sweeps A'awoop that each repUca performs. 



equilibrium value from opposite sides. This is illustrated 
in Fig. which shows qi and the quantity on the right- 
hand side (RHS) of Eq. (|13|l for an increasing number 
of sweeps for L — 128, T = 0.1. The two quantities do 
indeed approach each other from opposite sides, and once 
they agree they do not change if the simulation is run for 
longer. We applied this test for several lattice sizes using 
the parameters in Tabled and found that the data were 
well equilibrated in each case. 

We also require that the acceptance ratios for the 
global moves which interchange the different tempera- 
tures in the parallel tempering scheme be greater than 0.4 
on average and roughly constant as a function of temper- 
ature. Tabled shows the number of samples A'samp and 
the number of Monte Carlo sweeps A^swcop performed by 
each replica for each lattice size. The lowest temperature 
used is T = 0.10, the highest T = 1.40, wefi above the 
mean- field critical temperature — 1. 

We also use parallel tempering as an optimization al- 
gorithm to determine the ground state in the T — Q cal- 
culations. This is described in Sec. lVIII 



V. FINITE-SIZE SCALING IN THE CRITICAL 
REGION 

In this section, we determine the critical temperature 
Tc and the critical exponents for the disordered chain 
with a = 0.75. In Fig. [7| we show rescaled data for the 
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TABLE I: Parameters of the finite-T simulations. A'^samp 
is the number of samples, i.e., sets of disorder realizations, 
Aswcop is the total number of sweeps simulated for each of the 
2Nt replicas for a single sample, and Nt is the number of 
temperatures used in the parallel tempering method. 
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FIG. 7: The Binder ratio g as a function of L^^^iT - Tc) 
for several system sizes. We see that the data scale well for 
1/u ^ 0.30 ± 0.03 and Tc = 0.62 ± 0.03. This is a little 
lower than the (approximate) estimate of Bhatt and Young 
(Ref. I35I) who find Tc « 0.68. The inset shows the unsealed 
data for g. One can see that the data for different sizes cross 
at T - 0.62. 



Binder ratio g as a function of temperature for differ- 
ent system sizes [see Eq. (|1U|| ]. The data scale well for 
Tc = 0.62 ± 0.03 and 1/v = 0.30 ± 0.03, and agree with 
previous results by Leuzzi^ The inset of the aforemen- 
tioned figure shows the unsealed data for g which cross 
at Tc ~ 0.62. The error bars are estimated by varying g 
slightly until the data do not collapse well. 

We use the same approach below for the spin-glass 
susceptibility, for which the data, rescaled according to 
Eq. ((TT]) , are presented in Fig. |S1 The data scale well 
for r] — 1.34 ± 0.03 and the previously mentioned val- 
ues of 1/:^ and Tc- These values can be compared with 
the results of Kothar et al.'S^ r] = 3 — 2a (= 1.5), and 
1/v = (1/3) — 4e to first order in e, where e = a — 2/3. 
The latter gives 1/^ = suggesting that truncating the 
expansion to first order in e is not valid for a = 3/4. 

Having studied the critical region of the model to es- 
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Data for the rescaled spin-glass susceptibility 
for diff'erent sizes as a function of L^^-'iT - Tc). 
The data scale well for 1/u = 0.30 ± 0.03, Tc = 0.62 ± 0.03, 
and r] = 1.34 ± 0.03. 



timate T^ we are now able to simulate the system at low 
temperatures (T « 0.16Tc) in order to probe the spin- 
glass phase and test if 6* = 6*'. 



VI. RESULTS AT LOW TEMPERATURES 

In this section we study the model with cr ~ 0.75, d — 
1, marked by a cross in Fig.|21 at temperatures well below 
Tc ~ 0.62, in order to get information about large-scale, 
low-energy excitations in the spin-glass state. 

Differences between the various models for the spin- 
glass state can be quantified by studying^ji2i^£i2Si^ Pil), 
the distribution of the spin-glass overlap q. The RSB pic- 
ture predicts a nontrivial distribution with a finite weight 
in the tail down to q — 0, independent of system size. 
On the contrary, the droplet picture predicts that P{q) 
should be trivial with only two peaks at ±gEA, where qEA 
is the Edwards- Anderson order parameter. For finite sys- 
tems there is also a tail down to g = but this vanishes 
in the thermodynamic limit likei*^*^ ~ with 6' — 9. 

Figures 1^ and [TUl show data for P{q) at temperatures 
0.10 and 0.23, respectively. There is clearly a peak for 
large q and a tail down to g = 0. At both temperatures 
one sees that the tail in the distribution is essentially in- 
dependent of system size. A more precise determination 
of the size dependence of P(0) is shown in Fig.llllwhcre. 
to improve statistics, we average over the q values with 
l^l < qo, with qo = 0.50. The expectedi*^!^ behavior in 
the droplet model is P(0) L~^' , with 9' = 9 where^ 
9 = d — a. The dashed line in Fig. 1111 has slope —0.25, 
the expected value for a = 0.75 according to the droplet 
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FIG. 9: Data for the overlap distribution P{q) at T = 0.10. 
Note that the vertical scale is logarithmic to better make vis- 
ible both the peak at large q and the tail down to q = 0. In 
this, as well as in Fig. 1101 we only display some of the data 
points as symbols, for clarity, but the lines connect all the 
data points. This accounts for the curvature in some of the 
lines between neighboring symbols. In this paper, all distribu- 
tions are normalized so that the area shown under the curve 
is unity. 
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FIG. 10: Same as for Fig.|ibut at T = 0.23. 



model. The size dependence is consistent with a constant 
P(0), which implies that the energy to create a large ex- 
citation does not increase with size, at least for the range 
of sizes studied here. 

We perform a two-parameter fit of the data in Fig.llllof 
the form aL"'^' to find for T = 0.10, 9' = -0.005 ±0.009, 
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FIG. 11: Log-log plot of P(|q| < qo)/(2go), the spin overlap 
distribution for small q, against L for go = 0.50. The data 
are observed to not depend significantly on the system size L. 
The dashed line has slope —0.25, the theoretical prediction of 
the droplet model (Ref. 3). Asymptotically, the data should 
be parallel to this line according to the droplet theory. 

for T = 0.16, B' = -0.020 ± 0.007, and for T = 0.23, 
B' = -0.016 ± 0.006. The fitting probabihties"*! Q for 
the different temperatures are 0.135, 0.186, and 0.08, re- 
spectively, which is reasonable. Fixing d' — 0, we obtain 
fitting probabilities of 0.214, 0.053, and 0.076, respec- 
tively. On the contrary, fixing B' — 0.25, the predicted 
value for the droplet picture, we find Q < 10^''^ for all 
three temperatures. Hence we are confident that, within 
error bars, 0' ~ for the range of sizes studied. 

Moored has argued that should be trivial in the 
thermodynamic limit (i.e., just two delta functions at 
q = -^qEA) to first order in an expansion about the "lower 
critical dimension" cr = 1. This seems to disagree with 
our numerics, although possibly even larger sizes might 
be necessary to see the asymptotic behavior. 



VII. DOMAIN- WALL ENERGIES AT ZERO 
TEMPERATURE 

Since there are analytic predictions, Eqs. ||HJ) and ||SJ), 
for B for the present model, it is useful to directly cal- 
culate B from zero-temperature domain- wall calculations 
for different values of a. 

We use parallel tempering here as an optimization al- 
gorithm. We find that if we take the lowest temperature 
Tmin to be 0.05, then the minimum-energy state found 
at this temperature is the ground state. To test whether 
the true ground state has been reached, two criteria were 
adopted: (i) the same minimum-energy state has to be 
reached from two replicas at Tmin for all samples, and (ii) 
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TABLE II: Parameters of the T = simulations. The table 
shows the total number of Monte Carlo steps used for each 
value of a and L. 
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this state has to be reached during the first 1/4 of the 
sweeps in both copies. These conditions were satisfied for 
the parameters used in the simulations which are shown 
in Table in 

Having found the ground state with periodic (P) 
boundary conditions, we then flip the boundary condi- 
tions to "antiperiodic" (AP) by the following prescrip- 
tion. Consider one of the nearest neighbor (nn) bonds, 
e.g., that between sites 1 and L in Fig.^ Then, for all i 
and J, change the sign of the bond between i and j if the 
shorter path between those sites goes through the chosen 
nn bond. The new ground state is calculated and note is 
taken of the difference in energy 

5E = Eap - Ep (14) 

for each sample. On average, there is no preference for 
AP or P and so [^E'Jav — 0. Hence we take the average 
of the absolute value 

AE^[\SE\U (15) 

as a measure of the characteristic domain- wall energy. 

Data for a range of sizes and values of a are shown in 
Fig. El Our program automatically scales the interac- 
tions to set T^^ = 1, so no modifications are needed for 
the case a < 0.5 since the fact that the scale factor c{a) 
[in Eqs. and I0J) tends to zero for L — + cx3 in this 
range, is taken care of automatically. 

A smaller range of sizes has been studied for larger a as 
parallel tempering is less efficient in finding the ground 
state in this limit. At first sight this seems surprising 
because the nearest-neighbor interactions dominate and 
so there is no frustration, apart possibly from that due to 
the boundary conditions. The spins become stuck when 
a frustrated nn bond is larger in magnitude than the 
neighboring bonds. Parallel tempering does not seem to 
help much in this situation, perhaps because the local 
minima do not have a hierarchical structure. It would be 
interesting to investigate in more detail the conditions 
under which parallel tempering is efficient. 
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FIG. 12: Data for the defect energy AE =[\Ep~ Eap \ ]av for 
different sizes and values of a. The solid lines connect the data 
for a = 0.75, the value of a studied at finite temperatures. 



TABLE III: Fits for AE ~ aL" for different values of a. 
Q represents the quality of fit. The last column in the table 
shows the range of values of L used in the fits (for the possible 
values see Tableinj. For 0.10 < ct < 0.25, we exclude the data 
for L = 16 from the fits due to finite-size effects for small 
values of a. 

a a 9 Q L Range 






10 


-0 


818 


± 





033 





295 


±0 


007 





495 


32 - 


256 





25 


-0 


928 


± 





032 





301 


±0 


007 





759 


32 - 


256 





50 


-1 


023 


± 





033 





282 


±0 


005 





007 


16 - 


256 





62 


-1 


007 


± 





022 





239 


±0 


005 





235 


16 - 


256 





75 


-0 


925 


± 





022 





173 


±0 


005 





852 


16 - 


256 





87 


-0 


741 


± 





021 





079 


±0 


005 





801 


16 - 


256 


1 


00 


-0 


577 


± 





022 


-0 


023 


±0 


005 





319 


16 - 


256 


1 


25 


-0 


225 


± 





028 


-0 


239 


±0 


007 





222 


16 - 


128 


1 


50 





137 


± 





030 


-0 


457 


±0 


Oil 





349 


16 - 


64 


1 


75 





414 


± 





041 


-0 


647 


±0 


012 





150 


16 - 


64 


2 


00 





600 


± 





042 


-0 


788 


±0 


012 





727 


16 - 


64 


2 


50 





619 


± 





044 


-0 


913 


±0 


013 





289 


16 - 


64 


3 


00 





499 


± 





045 


-0 


935 


±0 


013 





975 


16 - 


64 



One expects that AE ~ L , where = 1 — cr in the 
long-range region 1/2 < cr < 2, and 9 = —1 in the short- 
range region cr > 2. Fits to the data obtained in Fig. 1121 
along with the analytic prediction, are shown in Fig. 1131 
and in Table IIIII In particular, there is a distinctive 
positive slope for cr = 0.75 in Fig. El (solid line), corre- 
sponding to a positive value of 6. This is in clear contrast 
to the finite-T data shown in Fig. El 

We see that the overall trends in the analytic prediction 
are reproduced by the data, but there is some round- 
ing where the analytical result has a change in slope. 
Note that our results for 6 always lie well below 1/2, the 
predicted^iSi value for cr = 1/2, and appear to saturate 
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FIG. 13: Variation of 9, with cr, obtained from fitting the 
data in Fig. ll2l to Ai5 ~ L". Also shown is the analytic pre- 
diction (dashed line): 9 = 1 — a for 1/2 < cr < 2, and 9 = — 1 
for (J > 2. The shaded region is where there is no thermody- 
namic limit unless the interactions are scaled appropriately 
by system size. 

at about 0.3 as cr —> 0, the SK limit. 

In recent work Aspelmeier et aZ.'*^. have given analytic 
arguments, based on the RSB theory, that the defect en- 
ergy should vary as A'''^ for short-range models in high 
dimension d, where N = L'^ and fj, is the exponent for 
the rms sample-to-sample fluctuations of the free energy 
of the SK model. A number of authorsMiiS4& have found 
that n ~ 1/4, at least at T = 0. We have therefore cal- 
culated for our one-dimensional model at T = for 
different values of a and show the results in Fig. In- 
deed we find that ^~l/4for(T^0 (the SK limit). Also 
note that it crosses over to the expected "self-averaging" 
value of II — 1/2 in the region cr > 1/2 where we no longer 
need to scale the interactions with system size. Since we 
work in d = 1 the expectation from Aspelmeier et alm^ 
is that 6 ^ for cr ^ 0. Figure [TH also shows 6 and, 
indeed, it is quite striking how 9 becomes close to /i at 
small cr. 



VIII. CONCLUSIONS 

We have performed finite-T Monte Carlo simulations 
and ground-state calculations of the long-range disor- 
dered Ising chain in one dimension with power-law in- 
teractions. For the ground-state calculations we have 
studied a range of values of a that encompasses all the 
possible regimes, from large a where the system is con- 
trolled by the short-range part of the interactions and 
has Tc — 0, all the way to near the SK model limit, 
cr = 0. For the finite-temperature Monte Carlo runs we 




a 



FIG. 14: Defect energy exponent 9 and the exponent for the 
rms sample-to-sample fluctuations in the ground-state energy, 
fj,, for different values of a. 



have concentrated on a single value of a, namely 0.75, 
which is in the range where there is a finite T^. 

For a = 0.75 we find the stiffness exponent for droplets, 
9', inferred from Monte Carlo simulations for sizes up 
to L = 512, to be close to zero, whereas the stiff- 
ness exponent for domain walls, 9, is positive. We find 
9 = 0.173 ± 0.005 (only statistical errors included) some- 
what lower than the analytic prediction^iSl that 6 = 0.25. 
These results are consistent with the RSB and TNT sce- 
narios according to which 9' = and 9^9'. Similar 
results have been found in earlier work on short-range 
models in three and four dimensions ^2*^^*^^ According 
to the droplet picture, 9' — 9 asymptotically, and there 
has been discussion as to whether the failure of the nu- 
merics to see this is an artifact of the small lattice sizes jii 
or whether it will persist in the thermodynamic limit. In 
the present work, the disagreement persists even though 
we are able to study a much larger range of L. 

Also noteworthy is our finding that 9 tends to a value 
of about 0.3 for cr ^ 0, the SK model limit. This is below 
the droplet theory prediction that 9 should tend to 1/2 
for cr ^ 0, but is close to the prediction of Aspelmeier et 
al^ that 6* ^ 1/4 in this limit. 
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